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Abstract 

We classify all integrable 3-dimensional scalar discrete quasilinear equations 
= on an elementary cubic cell of the lattice 1? . An equation Q3 = is called 
^ ■ integrable if it may be consistently imposed on all 3-dimensional elementary faces 



of the lattice Z'^. Under the natural requirement of invariance of the equation 
under the action of the complete group of symmetries of the cube we prove that 
the only nontrivial (non-linearizable) integrable equation from this class is the 
well-known dBKP-system. 
MSG: 37K10, 52C99 

Keywords: integrable systems, discrete equations, large polynomial systems, 
computer algebra, Reduce, Form, Crack. 



* Supported by the DFG Research Unit 565 "Polyhedral Surfaces" (TU-Berlin) 

^On leave from: Krasnoyarsk State Pedagogical University, Russia. SPT acknowledges partial 

financial support from the grant of Siberian Federal University (NM-project A^" 45.2007) and the 

RFBR grant 06-01-00814. 



1 



1 Introduction 



Although definitively shaped a decade ago, discrete differential geometry (see e.g. [3],ll]) 
has already provided much insight into structures that are fundamental both to classical 
differential geometry and to the theory of integrable PDEs. In addition to such purely 
mathematical fields, results in discrete differential geometry have a great potential in 
computer graphics and architectural design: it turns out that discrete surfaces param- 
eterized by discrete conjugate lines and discrete curvature lines — the basic structures 
in discrete differential geometry — have superior approximation properties and other 
useful features (see [9]). 
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Figure 1: lattice 

In this paper we consider the cubic lattice with vertices at integer points in 
the n-dimensional space = {{xi, . . . , € K}}. With each vertex (with integer 

coordinates {ii, . . . ,in), is £ ^) we associate a scalar field variable fi^,,,i„ G C. 

In what follows we need to consider the elementary cubic cell Kn = {{ii, . . . ,in) \ 
is G {0, 1}} of the lattice Z". The field variables /ji..i„ are associated to its 2" vertices. 
We will use the short notation f for the set (/00...0, • • • , /11...1) of all these 2" variables. 

An n-dimensional discrete system of the type considered here is given by an equation 
of the form 

Qn(f) = 0, (1) 

on the field variables on the elementary cubic cell Kn- For the other elementary cubic 
cells of the equation is the same, after shifting the indices of f suitably (see Figured]). 

In the last two decades the study of special classes of ([T]) which are "integrable" (in 
one sense or another) has become very popular. We give below only a brief account of 
the current state of this field of research, for a more detailed account cf. and the 

references given therein. In fact, discrete integrable systems underlie many classical 
integrable nonlinear PDEs, like the Krichever-Novikov equation and other examples. 
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the latter appear as a continuous limit along some of the discrete directions. Other 
well known classes of integrable geometric objects (with n = 3), like minimal surfaces, 
conjugate nets, constant curvature surfaces, Moutard nets, isothermic surfaces, orthog- 
onal curvilinear coordinates etc., are also obtained as some smooth limits along any 
two of the three directions of the respective discrete system. The remaining third dis- 
crete direction automatically provides us with a transformation known in the classical 
continuous geometric context as Jonas/Ribaucour/Backlund transformation between 
surfaces of the given class (see for more details). On the other hand, starting 

from the classical theorems on non-linear superposition principles and permutability 
of the aforementioned transformations between smooth surfaces of one of these types 
we obtain precisely the underlying discrete system. One of the cornerstones of the 
discrete differential geometry (the idea to look for cubic nonlinear superposition for- 
mulas of Backlund transformations of nonlinear integrable PDEs) was laid down in [3] . 
The duality between the smooth objects in any of the geometric classes of integrable 
smooth surfaces mentioned above and their Backlund-type transformations is therefore 
put into a symmetric form of a single discrete n-dimensional system and is encoded as 
the notion of {n + 1) -dimensional consistency [3]: 

An n-dimensional discrete equation (Cp is called consistent, if it may be imposed in 
a consistent way on all n-dimensional faces of a {n + 1) -dimensional cube. 

This can be also understood as the possibility to take and prescribe the n- 
dimensional equation ([1]) to hold on every ra-dimensional face of every elementary [n -\- 
l)-dimensional cube (of size 1, with edges parallel to the coordinate axes) without side 
relations to appear. For this reason ([T]) is often called a "face formula". A precise 
definition of consistency, suitable for the class of discrete equations treated in this 
paper, will be formulated in the next section. 

This paper is devoted to application of computer algebra systems Reduce and 
Form ([IO]), in particular the Reduce package Crack ([HI [12]), to the classification 
of 3-dimensional integrable discrete systems. 

The paper is organized as follows. In Section [2] we give a brief description of the 
known results on 2-dimensional integrable scalar discrete equations of type ([1]) and 
the precise definition of {n + l)-dimensional consistency condition for such discrete 
n-dimensional systems. 

Section [3] is devoted to the classification of symmetry types of quasilinear equations 
([T]) for dimensions n = 2, 3, 4. 

In Section m we describe the results of our computations (Theorem [2]): the only 
nontrivial (non-linearizable) integrable scalar quasilinear 3(i-face equation invariant 
w.r.t. the complete group of symmetries of the cube is given by the formula ([5]) below. 

Appendices A-E describe the technical details of the computations. 
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2 The setup 



The simplest but very important class of 2-dimensional integrable face formulas was 
investigated in detail in [H |2]. They have the form 

Q(/oo,/io,/oi,/ii) = 0, (2) 

where fij are scalar fields attached to the vertices of a square (see Fig. [2]) with two 
main requirements: 

1) Quasilinearity. ([2]) is affine linear w.r.t. every fij, i.e. Q has degree 1 in any 
of its four variables: Q = ci(/io, /oi, /ii)/oo + C2(/io, /oi, /ii) = C3(/oo, /oi, /ii)/io + 
C4(/oO) /oi) /ll) = • • • = ^iiii/oo/io/oi/ii + 9iiio/oo/io/oi + 91101 /oo/io/ii + • • • + Q'oooo- 

2) Symmetry. Equation ([2]) should be invariant w.r.t. the symmetry group of the 
square or its suitably chosen subgroup. 



/oil /ill 




Figure 2: Square K2. Figure 3: Cube K^. 



A few other requirements were given in [T], in particular the formula ([2]) involved 
parameters attached to the edges of the square. 

The second requirement of symmetry is obviously very important for the formu- 
lation of the condition of 3-dimensional consistency of ([2]). Namely, suppose we have 
an elementary cube (Sd-cell of Z^, cf. Fig. [3]) and impose ([2]) to hold on three "initial 
2(i-faces" {xi = 0}: Q(/ooo, /oio, /ooi, /on) = 0; {x2 = 0}: (5(/ooo, /loo, /ooi, /loi) = 0; 
{x3 = 0}: (5(/ooo, /loo, /oio, /no) = (these are used to find /on, /loi, /no from 
/0005 /looj fow, /ooi)- Then we impose ([2]) to hold on the other three "final 2(i-faces" 
{xi = 1}: (5(/ioo,/iio,/ioi,/iii) = 0; {x2 = 1}: Q(/oio, /no, /on, /in) = 0; {xs = 1}: 
Q(/ooi, /loi, /on, /ill) = 0; so for the last field variable /m we can find 3 (apriori) 
different rational expressions in terms of the "initial data" /ooo, /loo, /oio, /ooi- The 
3d- consistency is the requirement that these three expressions of /m in terms of the 
initial data should be identically equal. The subtle point of this process consists in the 
non-uniqueness of the mappings of a given square (Fig. [2]) onto the six 2(i-faces of the 
cube. The requirement of symmetry given above guarantees that we can choose any 
identification of the vertices of the 2-dimensional faces of the 3-dimensional elementary 
cube (Fig. [3]) with the vertices of the "standard" square where ([2]) is given; certainly this 
identification should preserve the combinatorial structure of the square (neighbouring 
vertices remain neighbouring). In [T] a complete classification of 3(i-consistent 2(i-face 
formulas (in a slightly different setting) was obtained; in [2] a similar classification was 
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given for the case when one does not assume that the formula ([2]) is the same on all 
the 6 faces of the 3-dimensional cube. 

In the next sections we give a symmetry classification of all possible Sd-face formulas 
defined on some "standard" 3-dimensional cube: 



with respect to the complete symmetry group of the cube. Here, as everywhere in 
the paper, indices of the field variables fijk give the coordinates of the corresponding 
vertices of the standard 3(i-cube where our formula ([3]) is defined. 

The requirement of consistency is now formulated similarly to the 2d-ca.se: given a 
4(i-cube with field values fijki, k, I G {0, 1}, one should impose the formula ([3]) on 
every 3d-ia.ce of it, by fixing one of the indices i,j, k, I, and making it for the faces 
which we will call below "initial faces", or respectively 1 for the faces which we will 
call "final faces" . One also needs to fix some mapping from the initial "standard" cube 
(with the vertices labelled fijk) onto every one of the eight 3(i-faces (for example {fnki} 
on {x2 = 1}). This can certainly be done using the trivial lexicographic correspondence 
of the type fijk fujk- Geometrically this lexicographic correspondence is less natural 
since it is not invariant w.r.t. the symmetry group of a 4d-cube. On the other hand 
there is an important example of such a non-symmetric formula corresponding to the 
discrete BKP equation (P, equation (76)). Another possibihty to avoid this problem 
is to impose the requirement of symmetry. More precisely, if one applies any one of the 
transformations from the group of symmetries of the Sd-cnhe, ([3]) shall be transformed 
into an equation with the left hand side proportional to the original expression Q: 
Q ^ \ ■ Q. Since this symmetry group is generated by reflections, one has = 1, 
so this proportionality multiplier A may be either (+1) or (—1) for any particular 
transformation in the complete symmetry group. 

From results in [T] we know that there are important 4(i-consistent "id-iace formulas 
which are preserved under a suitable subgroup of the complete symmetry group of the 
3(i-cube. No classifiaction of "id-iace formulas with such restricted symmetry property 
has beet caried out yet. 

3 Symmetry classification 

Every n-dimensional face formula Qn = which satisfies the requirement of quasilin- 
earity has a left hand side of the form 



with constant coefficients gu, where the summation is taken over all 2^" many 2"- 
tuples T) = (D00...O5 • • • 5 -D11...1), each power Di-^,__i^ of the respective vertex variable 
/n...i„ being either or 1. In other words: the 2" indices of are the exponents of 
the 2" vertex field variables each exponent -Dji...i„ being or 1. For example. 



QifoOO, /1OO5 /0IO7 /0OI5 /llO; /1015 /oil) /ill) — 



(3) 




(4) 
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Q2 = giiii/oo/io/oi/ii + giiio/oo/io/oi + giioi/oo/io/ii + • • • + ^0000 has 2^ = 16 terms, 
Qs has respectively 2^ = 256 terms, and (^4 has aheady 2^ = 65536 terms. 

In this Section we classify ra- dimensional quasilinear equations Qn = for n = 2, 3, 4 
that are invariant w.r.t. the complete symmetry group of the respective n- dimensional 
cube. This problem can be reduced to the enumeration of irreducible representations 
of this group in the space of polynomials of the form (j4]). Here this is done in a 
straightforward way: the group in question is generated by one reflection w.r.t the 
plane xi = 1/2 and [n — 1) diagonal reflections w.r.t. the planes xi = Xg, s = 2, . . . ,n 
(here Xk denote the coordinates in M"). 

To every reflection R from this generating set we assign (— ) or (+) and require 
the equality (5(f) = —Q{R{f)) (respectively (5(f) = Q{R{f))) to hold identically in 
all vertex variables f; this gives us a set of equations for the coefficients q-j). Running 
through all possible choices of the signs for the generating reflections we solve the 
united sets of simple linear equations for the coefficients for every such choice. The 
main problem consists in the size of the resulting set of equations: for n = 3 we have 
for each combination of ± for the 3 generating reflections around 770 equation for 
the 256 coefficients qDi...D8] for n = 4 every set of equations for the coefficients of Q4 
has around 250,000 equations for the 65536 coefficients q'Di...Di6- Naturally, not every 
combination of signs for the generating reflections is possible, most of the resulting sets 
of equations for gz)i...D(2n) allow only trivial solution gz)i...D(2n) = 0. The results of our 
computation are given in table form in Theorem [T] belowJl] 

As it turns out, for n = 2 three symmetry types of quasilinear expressions Q2 are 
possible. In the notation of Table [1] the flrst sign refers to the reflection on the line 
Xi = 1/2, the next sign stands for the reflection on the line xi = X2- For example, 
the expressions of the flrst type (H — ) are invariant w.r.t. to the reflection on the line 
Xi = 1/2, and show a change of sign after the reflection on the line Xi = X2- The last 
case (++) consists of expressions which are invariant w.r.t. any element of the complete 
group of symmetries of the square which corresponds to the choice of the (+) signs for 
the two generating reflections of the square w.r.t. the lines xi = 1/2 and xi = X2- 

For n = 3 alongside with the completely symmetric quasilinear expressions Q3 (the 

third case (+ + +) below), there are two other cases ( ) and ( — h +) of nontrivial 

quasilinear Q-^. In this notation the flrst sign refers to the reflection on the plane 
Xi = 1/2, the next signs stand for reflections on the planes Xi = X2, Xi = X3 (and 
Xi = X4 for n = 4). 

The resulting numbers of free coefficients gu in the symmetric face formulas Qn are 
given for each of the nontrivial cases in Table [H We also give the number of nonzero 
terms in Qn for each case. 

Especially remarkable is the totally skew-symmetric case ( ) for n = 3: it has 

^To solve the sparse but rather extensive hnear systems for the coefBcients appearing af- 
ter the sphtting w.r.t. the variables /ii...i„, a special linear equation solver had to be written. 
It can be downloaded together with other material related to this publication (for details see 
http: //lie .math.brocku. ca/twolf /papers/TsWo2007/readme|. 
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only one (up to a constant multiple) nontrivial expression 

Q( ) = (/loo /ooi)(/oio ~ /iii)(/ioi /iio)(/oii ~ /ooi)~ 

(5) 

(/ooi — /oio)(/iii — /ioo)(/ooo — /ioi)(/iio — /oil)- 

Precisely this expression gives the so called discrete Schwarzian bi-Kadomtsev-Petviashvili 
system (dBKP-system) — an integrable discrete system found in [6l [7] and studied in 
[T] where the fact of its 4(i-consistency was first established. The dBKP-system has 
many equivalent forms and appears in very different contexts. In addition to the 
known geometric interpretations and a reformulation as Yang-Baiter system ([5]), the 
dBKP-system may be considered as a nonlinear superposition principle for the classical 
2-dimensional Moutard transformations ([5]). 

The expression ([5]) enjoys an extra symmetry property: the equation Q( ) = is 

invariant under the action of the S'L2(C) group of fractional- linear transformations 

f ^ (af + 6)/(cf + rf) (6) 

(all group parameters a, 6, c, d are the same for all the vertices of the cube). This is 
in fact a direct consequence of its uniqueness in this class. Since this 5*^2 (C) action 
obviously preserves the symmetry type of an expression w.r.t. the group action of the 
cube symmetry group, it is reasonable to find the subclasses of S'L2-invariant face 
equations in each symmetry class. This is given in the third column of the table. 

Theorem 1 The nonempty symmetry classes of formulas Qn for face dimensions n = 
2,3,4 are: 



n 


types of symmetry, 

number of parameters and terms 


number of parameters 
in S'L2-invariant subcases 


2 


1) (H — ): 1 param.; 4 terms 

2) (— +): 3 param.; 10 terms 

3) (++): 6 param.; 16 terms 


1) 1 param.; 4 terms 

2) none 

3) 1 param.; 6 terms 


3 


1) ( ): 1 param.; 24 terms 

2) (—++): 13 param.; 186 terms 

3) (+++): 22 param.; 256 terms 


1) 1 param.; 24 terms 

2) none 

3) 3 param.; 114 terms 


4 


1) ( ) 

2) (+ ) 

3) (-+++) 

4) (++++) 


94 param.; 29208 terms 
77 param.; 26112 terms 
349 param.; 60666 terms 
402 param.; 2^^ terms 


1) 5 param.; 15480 terms 

2) none 

3) 3 param.; 15809 terms 

4) 18 param.; 96314 terms 



Table 1: Symmetry classification of the face formulas w.r.t the complete symmetry 
group of the cube 



7 



In particular, the explicit expression for the (non-S'L2-symmetric) 2(i-face formulas 

are: 



?i(/ii/io ~ /ii/oi ~ /lo/oo + /oi/oo) — Qiifu — /oo)(/io ~ /oi)) (7) 

- /lo - /oi + /oo) + g2(/ii/oo - /lo/oi) + (8) 
Qsifufwfoi ~ /ii/io/oo ~ /u/oi/oo + /lo/oi/oo); 

gi + g2(/ii + /lo + /oi + /oo) + qsifnfoo + /lo/oi) + (9) 

?4(/ll/l0 + /ll/oi + /lo/oo + /oi/oo) + 

?5(/ii/io/oi + /ii/io/oo + /ii/oi/oo + /lo/oi/oo) + 
Q'e/ii/io/oi/oo- 

The explicit form of the S'L2-symmetric face formula for n = 2 in the first (H — ) case 
is the same: 

Q = ?i(/ii - /oo)(/io - /oi)- 
For n = 2 and the symmetric case (++) the 5'L2-symmetric face formula is: 

Q = Qiifnfw + /ii/oi — 2/11/00 — 2/10/01 + /lo/oo + /oi/oo)- 

It should be noted that ([8]) is 3(i-consistent; it can be simplified using SL2 trans- 
formations and put into explicitly linear form Q = fn — fio — foi + /oo or into 
Q = /ii/oo — /10/015 so the respective equation Q = is equivalent to a linear equation 
log /ii + log /oo = log /lo + log foi- 



(+-): Q = 
(-+): Q = 

(++): Q = 



4 Nonexistence of nontrivial integrable face formu- 
las in other symmetry classes for n = 3 

In this section we give a computational proof of our main result: 

Theorem 2 Among the three possible symmetry types of 3- dimensional quasilinear 
face formulas given in the second column of Table [21 only formula ^ gives a non- 
trivial 4:d- compatible face equation. Any Ad- compatible face formula in the other two 
symmetry types may be transformed using the action of the group 5/72 (C) on the field 
variables to one of the following linearizable forms: 

Q^^^ = /ooo/ooi/oio/oii/ioo/ioi/uo/iii ~ (10) 

Q^'^^ = /oOl/oio/lOo/lll ~ Cr/ooo/oil/lOl/llO) (11) 

Q^^^ = (/ool + /olo + /loo + /ill) - ^(/ooo + /on + /loi + /no), (12) 
where a = ±1. 

Technically, in order to find 4(i-consistent 3(i-face formulas among the other two 3- 
dimensional cases ( — 1-+) and (+++) (as listed in Tabled]), one shall run the following 
algorithmic steps: 
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Step 1. Take a copy of this Q3 formula, map it onto the four initial faces of the Ad- 
cube (where one of the coordinates Xj = 0), solve the mapped equations with respect 
to /oiii, /loii! /iioi and /mo leaving the other variables free. 

Step 2. Then substitute the obtained rational expressions for /om, /loii, /iioi and 
/mo into the copies of the face formula mapped onto the four final faces (where one 
of the coordinates Xj = 1), finding respectively four different expressions for the last 
vertex field /mi. 

Step 3. Equate these 4 expressions for /mi obtaining three rational equations in 
terms of 11 free variables /oooo, /ooio, /oioi, • • • and the parametric coefficients qDi...Ds 
left free in the chosen symmetry class. 

Step 4- Removing the common denominators of the equations and splitting the 
resulting polynomials w.r.t. the 11 free variables fijki one obtains a polynomial system 
of equations for the free coefficients qDi...Ds- 

Step 5. The latter should be solved, resulting in a complete classification of Ad- 
consistent quasilinear scalar 3(i-face formulas. 

This approach, applied in a straightforward way, results in extremely huge expres- 
sions. Even building the rational expressions in Step 4 in a straightforward way seems 
to be unrealistic: for a typical Sd-face formula from Table [1] Step 4 should end up 
(as our test runs allowed for an estimate) in an expression with around 10^^ terms, 
which is beyond the reach of computer algebra systems in the foreseeable future. Even 
brute force testing of 4(i-consistency of the smallest solution (face formula ([5]) which 
has no free parametric coefficients qDi...Ds) results in ^ 2 ■ 10^ terms (after substitut- 
ing the expressions for /om, fion, /iioi; /iiiO; collecting the terms over the common 
denominator and expanding the brackets before the cancellation can start in Step 
Technically this is explained by the presence of 4 different symbolic denominators of 
the rational expressions for /om, fwu, /iioi, /iiio and their various products. A care- 
ful step-by-step substitution and cancellation of like terms in several stages still can be 

done even on a modest computer for this ( ) case. Using the system Form (this 

system was specially designed for large symbolic computations), one can prove that all 
terms finally cancel out for the case of the integrable 3d-ia.ce formula ([5]) thus giving 
a computational proof of its 4(i-consistency in 3 min CPU time (3 GHz Intel running 
Linux SUSE 9.3) and less than 200 Mb disk space for temporary data storagej^l 

As the estimates given in Appendix A show, the straightforward approach based 
on Steps 1-4 is unrealistic for the other two 3-dimensional cases listed in Table [U even 
at the stage of generation of the consistency conditions (Step 4). 

In order to classify discrete integrable 3(i-face formulas Q3 = for the case ( — h-l-) 
and the hardest case (-I--I--I-) we used a totally different randomized "probing" strategy, 
explained in detail in Appendices B, C, D. 

After the computation (cf. Appendices B-E) the list of candidate face formulas Q3 
for the case (+++) included 5 formulas (before the verification that these formulas. 



^The commented Form logfile for this estimate can be downfoaded from 



http: //lie .math.brocku. ca/twolf /paper s/TsWo2007/BruteForceCh.eck/Case- 

•^The Form code of this run and its logfile can be downloaded from 



http : //lie .math.brocku. ca/twolf /papers/TsWo2007/BruteForceCheck/Case /). 
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obtained by our "probing" method, really give 4(i-consistent 3(i-formulas). For the case 
( — the list of candidate face formulas included 3 formulas. All of them include a 
few free parameters. As one can show, all these formulas can be greatly simplified using 
the action of the group SL2{C) on the field variables f, resulting in the 4(i-consistent 
3(i-face formulas (ITOl) . (ITTl) . (fT2l) . The first two formulas (ITOl) . (fTTl) can be linearized 
using the logarithmic substitution fij^ = log/j^fc. 

The expressions for the aforementioned candidates, the Form procedures and their 
logfiles showing the simplification process can be downloaded from 
[http : / /lie .math . brocku . ca/twolf /paper s/TsWo2007/ . Here we just give one ex- 
ample of such a 5'L2-simplification: the Q-expression for the case Ql in 



http : //lie .math. brocku. ca/twolf /paper s/TsWo2007/SL2-simplifi cat ion/Case+++/ 



IS 



Q = 9io5(/ooi/oio/ioo/iii + /ooo/oii/ioi/iio) + 

?io7(/ooi/oio/ioo/iio/iii + /ool/oio/ioo/ioi/iii + /ool/oio/oii/ioo/iii + 
/ooo/oii/ioi/iio/iii + /ooo/oii/ioo/ioi/iio + /ooo/oio/oii/ioi/iio + 
/ooo/ooi/oii/ioi/iio + /ooo/ooi/oio/ioo/iii) + 

^^(/ooi/oio/ioo/ioi/iio/iii + /ool/oio/oii/ioo/iio/iii + 

/ool/oio/oii/ioo/ioi/iii + /ooo/oii/ioo/ioi/iio/iii + 
/ooo/oio/oii/ioi/iio/iii + /ooo/oio/oii/ioo/ioi/iio + 
/ooo/ooi/oii/ioi/iio/iii + /ooo/ooi/oii/ioo/ioi/iio + 
/ooo/ooi/oio/ioo/iio/iii + /ooo/ooi/oio/ioo/ioi/iii + 
/ooo/ooi/oio/oii/ioi/iio + /ooo/ooi/oio/oii/ioo/iii) + 

-2~(/ooi/oio/oii/ioo/ioi/iio/iii + /ooo/oio/oii/ioo/ioi/iio/iii + 

?105 

/ooo/ooi/oii/ioo/ioi/iio/iii + /ooo/ooi/oio/ioo/ioi/iio/iii + 
/ooo/ooi/oio/oii/ioi/iio/iii + /ooo/ooi/oio/oii/ioo/iio/iii + 
/ooo/ooi/oio/oii/ioo/ioi/iii + /ooo/ooi/oio/oii/ioo/ioi/iio) + 

2^^(/ooo/ooi/oio/oii/ioo/ioi/iio/iii) 

9l05 

The simplifying transformation consists in the following steps: 

4 

1- fijk ^ ^^fijk and Q ^ (this ehminates the parametric gios and gio?) 

2. fijk ^ -J— (and removing the denominator in Q afterwards) 

3- fijk ' fijk 1- 

This produces the simplified form Q = /ooi/oio/ioo/iii + /ooo/oii/ioi/iio- 
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Appendix A: The Size of Consistency Conditions 

The following considerations are made under the assumption of generic unknown coef- 
ficients q-x) in the face formula (HI) not satisfying additional symmetry conditions. 

Any {n + l)-dimensional hypercube built from 2"+^ vertices fij^..A„+i, ik £ {0, 1} has 
2(n + 1) faces located in the (logical) planes = and Xk = 1, k = 1, . . . ,{n + 1). 
The face relations for the n + 1 faces that correspond to = are 

= E^^ n (/n.....-.o.,,,....„,0'''-"=-^"=- -^-^ (13) 

D is=0,l 

They can be used to determine /i,,ioi..i with the being in the k^^ index position. Each 
of these face relations involves 2" /-variables and thus 2^" terms, half of them include 
/1..101..1 as a factor and the other half not. Solving the face relation Xk = 0: 

= Afc/i,.ioi..i + Bk (14) 

where A^, are expressions in gx), fp for /i,,ioi..i and substituting /1..101..1 = —B^/Ak 
in any expression that involves /1..101..1 linearly (like other face relations) and taking the 
numerator over the common denominator amounts to multiplying all terms that involve 
/1..101..1 by —Bk and all other terms by Ak- As Ak and Bk involve each 2^"/2 = 2^"~^ 
terms this means that a substitution of /1..101..1 increases the number of terms by a 
factor of 2^"~^, before cancellations and reductions will be made. 

The 2°*^ half of face relations for the n + 1 faces that correspond to = 1 are 

= E^^ n (/n....-.n......„.j''" (15) 

D is=0,l 

Each one of them involves and n of those /-variables which have exactly one as 
index in any one of the n + 1 index positions apart from the k^^ position. Replacing each 
one of these n /-variables by using the corresponding = face relation increases the 
number of terms by a factor 2^""^ each time, giving in total 2^" (2^"~^) = 22"("+^)~"' 
terms. In each substitution the degree of the coefficients qd increases by one, reaching 
finally n + 1. 

Solving one of the n + 1 many x^ = 1 face relations 

= Gkfu.A + Hk (16) 
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dimension of face 


n 


2 


3 


4 


5 


# of /-variables in face formula 




4 


8 


16 


32 


# of terms in face formula 
(= # of undetermined 
coefficients in Qn ) 


22" 


16 


256 


65536 


4.3 X 10^ 


# of all /-variables 

in (n -I- l)-dim. hypercube 




8 


16 


32 


64 


# of indep. /-variables 
in (n -I- l)-dim. hypercube 


2"+i - n - 2 


4 


11 


26 


57 


# of n-dim. faces 

in (n -I- l)-dim. hypercube 


2(n + 1) 


6 


8 


10 


12 


# of consistency conditions 


n 


2 


3 


4 


5 


upper bound on the # of terms 
of each condition 


2{2" + l(ii+l)-2n-l} 


5.2 X 10^ 


1.4 X 10^^ 


2.8 X 10^^ 


1.9 X 10"2 


total degree of the 
in each condition 


2n -1- 2 


6 


8 


10 


12 


upper bound estimate of the 
# of equations resulting from 
splitting each condition 


2n(2n + 2)(2"+'-"-3) 


864 


6.4 X 10^ 


8.0 X 10^5 


2.7 X 10^1 


estimated average # of terms 
in each equation 


22"+l(n + l)-2n-l 
2n(2n+2)(^"^^"""^) 


606 


2.2 X 10^ 


3.5 X lO^'' 


7 X 10-'^'' 



Table 2: Some statistics of faces and consistency conditions 



for /11..1 and substituting /n 1 = —Hk/Gk in the other face relations gives n indepen- 
dent consistency conditions 

GjHk = GkHj, J = l,...,k-l,k + l,...,n + l (17) 

with each Gi and Hi having 22"("+^)~"/2 terms, i.e. each consistency condition involving 
2^2'^"-{n+i)-n-iY terms. The total number of terms of the n consistency conditions is 
thus n2{2"^^("+i)-2"-i}. 

To compute an upper bound of the number of conditions that result from splitting 
each consistency condition with respect to the independent /-variables we note that 
their highest degree is equal to the total degree of all q-ri, i.e. it is 2n + 2. The only 
exception is /00...0 which does not occur in the = Qn\i^k face relations and enters 
only through substitutions, so its highest degree in the constraints is 2n. We thus get 
for an upper bound of the number of different products of different powers of /00...0 and 
the other 2"+^ — n — 3 independent /-variables the value 2n(2n + 2) ^ -n-3) ^ -^j^j^ 
this number and the number of terms of each constraint we get with their quotient an 
estimate of the average number of terms in each equation (see Table 2). 

Remark. Although, strictly speaking, these are upper bounds for the size of condi- 
tions, one shall keep in mind that any computer algebra system shall generate all these 
terms expanding the brackets and only after generating all of them or a part of them 
it can search for possible cancellations or reductions. As our test runs with Form 
had shown, for the 3-dimensional case (+ + +) given in Table [1] (with a much smaller 
number of independent qd than 2^" but the same number of terms in Q3), the total 
number of terms to be generated for each of the 3 consistency conditions is around 10^^ 
(compared to 1.4 ■ 10^^ in Table 2). A typical single equation for the coefficients qxi 
resulting from splitting a partially formed consistency condition has a few thousand 
terms of degree 8 in the parametric gD- 
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Appendix B: The computational problem 



As outlined before (cf. Section Hj), the task consists of the following steps. 

1. Formulate the relations for the n + 1 faces Xk = 0. 

2. Solve them for /ii...ioi...i- 

3. Formulate the relations for the n + 1 faces = 1. 

4. Perform the substitutions obtained under 2. in the relations of 3. 

5. Solve one of the resulting relations for and 

6. Substitute fu,,,i in all other n face relations of 4. 

7. Split these consistency conditions with respect to all the occurring independent 
/-variables to obtain an overdetermined system of equations for the unknowns 
q-D- 

8. Find the general solution of this system. 

9. Reduce the number of free parameters of the solutions using 5'L2(C)-transformations 
©. 

Although the algebraic system for the unknown coefficients is heavily overdeter- 
mined the following difficulties appear. 

1. Strictly speaking, in order to formulate even only the smallest subset of conditions 
one would have to formulate at least one consistency condition (by performing 
steps 1., 2. fully and 3.- 6. for at least two Xk = I face relations f|T5|) before step 
7) i.e. to generate an expression with 2'^^"^ (n+i)-2n-i} i^gj^j-Qg^ 

2. If one found a way around this hurdle then the resulting equations are of high 
degree 2n -|- 2 with on average many terms. 

3. Even if one were able to generate 100,000's of equations and thus find shorter 
ones which one could solve for some unknowns in terms of others, one would face 
the problem that many cases and sub-sub-cases have to be investigated due to 
the high degree of the equations and 

4. that one has to generate billions of equations to find some that are independent 
of the ones generated so far. As we explain below, one can not hope that the 
first, say, 10^ conditions will be equivalent to the full system of equations for qxi, 
even though we have only very few unknown g^'s due to the "triangular" form 
of this huge system as we explain in Appendix D. 

In the computations to be described in this paper we have n = 3. The full problem 
(without cubical symmetry) would mean to generate 3 consistency conditions involving 
each about 10^^ terms that split into an estimate of 10^° polynomial equations each 
being homogeneous of degree 8 for 256 unknowns and involving on average over 10^ 
terms. 

To make progress we introduce different solving techniques but also simplify the 
problem: 
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1. We restrict our problem to face formulas that obey the full cubical symmetry (cf. 
Section [3]). This reduces the problem to 3 cases: 

• case ( ) with 232 identically vanishing and 24 other q-j) depending on 

1 parameter, 

• case ( — h +) with 70 identically vanishing and 186 other q-x) depending on 
13 parameters, and 

• case (+ + +) with no identically vanishing gx) and all 256 gx) depending on 22 
parameters. 

In the third case (+ + +)- the hardest case - the fact that only 22 of the gxi are 
parametric simplifies the problem but the simplification is limited because none 
of the other 234 g^i needs to vanish just because of the extra cubical symmetry 
(see Section [3]). 

2. We ease the formulation of the large consistency conditions (before splitting them 
into smaller equations) by replacing z of the v := 2"+^ — n — 2 independent /- 
variables by zero, and replacing u of them by random integer values wn-equal 
zero, leaving s := v — z — n of them in symbolic form. As a consequence, we 
only get a set of necessary and not sufficient equations for the g^ but we can 
repeat this procedure on a gradually increasing level of generality (by increasing 
s and lowering z). A crucial feature is that the preliminary knowledge about the 
solution is used in the formulation of new necessary systems of equations (see 
Appendix C). 

3. We design a dynamic process that automatically organizes an iteration process 
that generates and solves/simplifies/'makes use of necessary equations automat- 
ically. As an important feature, a set of newly generated necessary equations is 
not read into the ongoing computation at once but gradually on demand (see 
Appendix D). 

4. A recent extension is the parallelization of different case investigations on a com- 
puter cluster. 

We discuss points 1-3 now in more detail. 

Appendix C: Random probing 

In this Appendix we explain the "probing technique" in more detail which was men- 
tioned briefly in Appendix B under point 2. in the list of techniques. It is based on 
replacing a number z of the v := 2""*"^ — n — 2 independent /-variables by zero (for 
n = 3, f = 11), and replacing a number u of them by random integer non-zero values, 
leaving s := v — z — u of them in symbolic form. 

The computation has 3 phases: finding solutions, verifying the obtained solutions 
probabilistically and again rigorously. The first two phases use the probing technique. 
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Still, it makes sense to distinguish them because optimal values of u and s differ in 
both cases. 

In the probing technique the two types of replacements, i.e. replacing by zero 
or by a non-zero integer, share the same disadvantage: the number of independent 
parameters, i.e. of symbolic fa which allow to split the consistency condition into 
many smaller equations, is reduced by one. 

The advantage of replacing an fa by zero instead of a non-zero integer is that 
expressions shrink more. Also, replacements by non-zero integers often introduce extra 
solutions to the generated conditions and it appears to be costly to eliminate these 
spurious temporary solutions by leading them to a contradiction with more conditions 
to be generated based on other random replacements. Therefore it is more productive 
to replace, for example, one fa by zero than to replace two fa by non-zero integers. 

Consequently, in the first phase of finding solutions at most one replacement by a 
non-zero integer is made (i.e. u < 1) and after starting with z = 9,u = 0, s = 2 one 
increases generality gradually by either changing u from to 1 and decreasing 2 by 1, 
or by decreasing u from 1 to and increasing s by 1 until z = l,u = 0, s = 10. We 
need the occasional substitution by one non-zero integer because we want to increase 
the generality in as small as possible steps in order to avoid the generation of too many 
too high degree equations with too many terms. This would happen if we decrease z 
by 1, keep u = and increase s by 1. A run in full generality z = 0,u = 0, s = 11 
is computationally prohibitive, therefore in a second phase of confirming the found 
solution probabilistically one starts with z = 4,u = 1, s = 6 and increases generality 
by decreasing z, increasing u and keeping s constant until z = 0,u = 5, s = 6. By 
testing a hypothetical solution with this final setting many times, the correctness of 
the solution is confirmed with an arbitrarily high probability. 

In both phases we want to generate as few and as simple equations as possible in 
each generation step. So we continue using the same setting of z, u, s as long as possible 
(i.e. as long as there are still resulting new conditions after randomly choosing other 
sets of z many to be 0, s of fa to be kept symbolic and randomly assigning integer 
non-zero values to the other u parametric fa) before generalizing it, i.e. making s larger 
and z smaller. This is regulated by one parameter which specifies the maximum number 
of consecutive times that a 'probing' (generation of conditions) attempt yielded only 
identities before changing z, u, s. 

In a third phase, after all hypothetical solutions have been obtained and been 
checked probabilistically, they are checked again, now rigorously. This has been done 
either by a brute force check using the computer algebra system Form or by using 
S'L2(C)-transformations on the field variables f to reduce solutions to integrable trivial 
forms (UnD, (mD, ([I2D in Section H 

A helpful and initially unexpected feature of the probing technique is that the 
resulting equations appear to be somehow triangularized in the following sense. Each 
unknown qiii2. rn = 2", ij G {0, 1} is the coefficient of a product of Zi + Z2 + • • • + 
many different factors /„. That means, that at the beginning of the computation when 
many fa are replaced by zero, the q-ri with a high index sum do not occur. Only later 
on as fewer fa are replaced by zero, gradually q-n with higher index sum appear in 
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the equations. On one hand this is a good feature, providing a partially triangularized 
system of equations. On the other hand this means that although we only want to 
compute a relatively small number of unknowns (for n = 3 and case (+++) these are 
22 q-j)) it is not enough to formulate only a comparable number of the huge total set of 
equations (about 6.4 x 10^ equations for n = 3). A set of equations that is equivalent 
to the complete set of equations is only obtained towards the end of generalizations. 
For example, as the coefficient of the product of all the 2"' many f occurring in 

at least one face formula, at most 2" — 1 of them belonging to the independent /„, 
can only occur in at least one consistency condition if none of the fa occurring in at 
least one face formula is replaced by zero, i.e. if u and s are big enough to satisfy 
{u + s)> (2"-l). 

That in turn means that millions of the early equations are redundant which implies 
a large inefficiency in generating equations. This can be avoided by using known 
relations q-j) = hxi^q'r^), which were derived in the solution process so far, as automatic 
simplification rules when generating new equations. 

Three problems remain to be considered. 

1. By replacing fa through numbers it may happen that Ak in any one of the 
Xfc = face relations (fT^ becomes zero. Then a new equation generation attempt with 
different or more general random replacements has to be made. 

2. Similarly, it may happen that the coefficient of qii,,,i in all + 1 many = 1 
face relations f|T6|) is zero. Then a different replacement has to be tried as well. 

3. Even if none of the in f|T6|) becomes zero it may happen that and -Bfc in 
one face relation are not prime and then a solution for the gn which makes the greatest 
common divisor GCD{Ak, B^) to zero is potentially lost when performing substitutions 
/1..101..1 = —Bk/Ak- The same applies to common factors of Gk and Hk in the single 
face relation (fT6|) that is applied to replace Therefore each consistency condition 
has to be multiplied with a product of all common factors of any pair A^, Bk and of 
all common factors of the pair Gk-, Hk which is used to replace To lower the 
computational cost one drops multiplicities of the factors. These factors involve in 
general as well as fa and therefore the multiplication has to be done before splitting 
the consistency condition with respect to the independent /„. Alternatively one can 
split the consistency condition before multiplication and instead multiply and duplicate 
the equations in the following way. 

Let P{qD, fa) be one of the above mentioned factors and let = P be split into 
a system = Pilqi)) where redundant equations are droppedQ Instead of multiplying 
a constraint = C{qxi,fa) with P, splitting the equation = PC into individual 
equations and factorizing all of them afterwards, it is equivalent but much more efficient 
to split = C into a system of equations = Cjiqi)) and = P into a system 
= Piiqx)) and to consider the equivalent system = PiCj, 'ii,]. 

To summarize, in order not to loose solutions for the gn, the procedure is 

^The definition of 'redundant' depends of the effort one wants to spend at this stage. In the 
implementation of this algorithm the polynomials Pi are divided by the coefficients of their leading 
terms with respect to some ordering of the and then duplicate Pi arc dropped. 
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• to collect all common factors Pr of all pairs Ak-, Bk and of the pair Gk-, Hk used 
to substitute 

• to drop duplicate factors, 

• to split all consistency conditions giving a system S of equations = Cj, and 

• to split for each factor Pj. the equation = Pr into a system = Pn-, i = 1 . . .i^ 
where again equations that are redundant within one such system are dropped. 

• If a system = Pri includes a non- vanishing Pri either because Pri = 1 or because 
PriiQTi) is known to be non- zero based on the inequalities that are known for some 
qxi then this system is ignored because the corresponding Pr is non-zero. For every 
other such system = Pri, replace the system S of conditions = Cj by the new 
system 5" consisting of the equations = PriCj, 

Appendix D: A Succession of Generating and Solving 
Equations 

The system of algebraic equations for the is investigated by the computer algebra 
package Crack that aims at solving polynomially algebraic or differential systems, 
typically systems that are overdetermined and very large. It offers various degrees of 
interactivity from fully automatic to fully interactive. The package consists of about 40 
modules which perform different steps, like substitutions, factorizations, shortenings, 
Grobner basis steps, integrations, separations,... which can be executed in any order. 
In automatic computations their application is governed by a priority list where highly 
beneficial, low cost and low risk (of exploding the size of equations) steps come first. 
Modules are tried from the beginning of this list to its end until an attempt is successful 
and then execution returns to the start of the list and modules are tried again in that 
order. This simple principle is refined in a number of ways. For more details see [TT] , 

In order to accommodate the dynamic generation of equations and their successive 
use, all that had to be done was to add two modules: 

• one for generating a new set of necessary conditions using the 'probing' technique 
from Appendix C and writing the generated equations into a buffer file, and 

• another module for reading one non-trivial equation from this buffer file (i.e. 
for continuously reading equations until one is obtained that is not instantly 
simplified to an identity modulo the known equations or until the end of this file 
is reached), 

and to determine the place of these modules in the priority list. The need for a buffer file 
arose because the number of equations generated in each 'probing' is unpredictable. 
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especially in view of the large impact that factors Pri can have on the number of 
equations (see the end of Appendix C)!^! 

Some more miscellaneous comments. 

As shown in Appendix C the generated equations often take the form of products 
set equal to zero. This leads to many case distinctions of factors being either zero 
or non-zero and consequently other factors being zero. The depth of sub-case levels 
sometimes reaches 20. Because buffer files are only valid for the case in which they were 
generated (because they make use of case- dependent known substitutions = hq^^q'r^) 
and case-dependent inequalities) and because of these deep levels of sub-cases, the 
number of buffer files easily reaches 100,000 and more (e.g. too many to be deleted 
with the simple UNIX command rm *). Therefore the case label is encoded in the 
buffer file name allowing to delete buffer files automatically when the case in which the 
file was created and all its sub-cases are solved. 

There is much room for experimenting with the place of the two new modules within 
the priority lisl§. On one hand one wants to read and create early, so that the ongoing 
computation has many equations to choose from when looking for the most suitable 
substitutions, shortenings, ... . The problem is that these equations are all generated 
with the same limited information on relations between g^, and thus they have a high 
redundancy. Also, dealing with many long equations does slow down Crack. On the 
other hand, giving the branching of the computation into sub-cases a higher priority 
generates an exponential growth of sub- and sub-sub-cases which drastically increases 
the number of buffer files to be generated because they are only valid for the case they 
were generated for or for its sub-cases. 

With each investigated case, say ^5 = 0, the other case 7^ generates inequalities 
which the package Crack collects, updates and makes heavily use of to avoid further 
case distinctions as far as possible, and in this computation also to drop factors of the 
Pmi as mentioned in Appendix D. 

The individual cases can either be investigated serially or in parallel. 

If two solutions of two different cases, for example, the solutions for ^5 = and 
gs 7^ can be merged into one analytic form, if necessary by a re-parametrization, 
then this is achieved by one of the modules of Crack ([13]). 

Appendix E: The Computation 

The computation was not performed in a single run. It fact due to the big workload 
it extended over several months. Simple subcases were solved initially whereas harder 

^This arrangement is similar to the design of CPU chips which have not only access to their own 
register memory the equations known within Crack) and access to the hard disk the possibility 
to call the 'probing' module) but which also have access to cache memory (~ the buffer file). 

^apart from the necessity to give the 'reading from the buffer' module a higher priority than the 
'creating of a buffer' module in order to try reading and emptying a buffer file first (if available and 
not already read completely) and to create a new buffer file only if none is available or if the available 
one is already completely read in 
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ones were completed only after the probing technique from Appendix C and its au- 
tomatic interplay with the package Crack were developed. Even then it took some 
time to fine-tune parameters and put the new modules in the right place within the 
priority list of procedures in Crack. Finally, it was nearly possible to do the com- 
putation fully automatically, only a few times the proper case distinctions had to be 
initiated manually at the right time to be able to complete the computation. If one 
would add up purely the necessary computation time without runs following a poor 
manual choice of case distinctions leading to computations which generated too large 
systems and which could not be completed then this would amount to 2 weeks of CPU 
time on a 3GHz AMD64 PC. 
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